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A6STEj!iG-I 


The present work describes a methodology for solving 
tne uncapacitated warehouse location problem considering 
demands as independent random variables belonging to a 
known 2-par8meter distribution. The objective is to 
minimize total fixed and variable transportation costs. The 
problem is niitially structored as a chance coist rained 
problem assuming that the demands are normally distribuxed. 
The chance constrained programjTung problem is converted 
into a deterministic mixed integer programming model. The 
mixed integer piogranming model is solved to determine if 
aJ-1 the decision variable are integers. Integer values for 
all the decisiOii variable indicate that an optimal solution 
to the problem has been obtained. If all the decision 
vaiiables are not integer the proble'i' is solved using branch 
and bound concepts. An efficient branch and bound algorithn 
involving three node simplification procedures and eight 
branching rules has been presented to arrive at the optimal 
solution. The branch and bound algorithm used is tested 
using a problem from the literature. The computational 
efficiency of the algorithm has been evaluated based on the 
comxjutational expeiience of solving 20 problems. In all. 



problems of four different sizes are considered pjid for 
eac-i size five problems r re solved. It is found that the 
computational tiiCe does not increase drastically with the 
increase in the size of the problem. 

A computer package for the proposed branch and b^und 
procedure has been developed for the I3!I 7044/1401 system. 

The package in its present form can handle problems involving 
warehouses and lOO customers. A procedure for the selection 
of optimum confidence level with which the consumer demand 
will be satisfied lias been developed and illustrated with 
the help of a numerical exa-Hi^le. 



CHAPTEE I 


lUTROEUaTIOE 


Management of industrial logistics seeks to maximize 
the economic value of products or materials hy getting and 
having the products where they are wanted, at the time they 
are wanted and at a reasonable cost. The management of 
physical supply and distribution forms an important subsystem 
of the overall domain of the management of industrial logis- 
tics and considerably influences the cost of the product/ 
commodity to the consumers. The design of a physical distri- 
bution system involves the selection of the number, size and 
location of warehouses, the mode of transportation from plant 
to the warehouses and from warehouses to the customers. It 
is interesting to point out that operation of warehouses and 
carrying cost of material stored in them account for an amount 
equivalent to 6 percent of gross national product. On the 
individual firm level the cost of physical distribution of 
the product may be as high as 15 percent of the total expen- 
diture, These figures reflect the importance of designing 
an efficient-distribution system. 
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During the past few years a large number of organi- 
zations, attracted by the possibility of making substantial 
savings, through the use of efficient distribution systems, 
have been critically examining in warehousing policy. 

Since the selection of location and size of the warehouses 
considerably influences the cost of distribution, a lot of 
attention is currently being devoted to this specific area. 

There are a host of tangible and intangible factors 
which influence the location of a warehouse. Till the 
emergence of operation research techniques, the problem 
lay mainly in the hands of economists who handled it 
primarily using qualitative approaches. However, there 
are certain quantitative aspects of the problem which cannot 
be reckoned accurately through intuition alone. Since 1958, 
considerable research effort seems to have been deployed 
for the development of quantitative approaches for warehouse 
location problems. Though, today quantitative techniques 
are heavily emphasised, one should not overlook the importance 
of qualitative aspects of this problem. It should be realized 
that analytic al approaches yield a solution to the model but 
not necessarily to the problem. The solutions obtained 
through analysis serve only as an aid in decision making. 
There remain a number of non-quantifiable factors which 
must be considered alongwith analytical results in making 



the transition from the model to the problem. In some cases 
the solution obtained from analysis will appear impractical 
from an operational point of view. Such solutions should 
be interpreted as bench marks against which operationally 
accepted solutions can be compared. 

The selection of the set of potential warehouse locations 
depends on many tangiDle and intangible factors. Important 
factors which influence the location of a warehouse are; 
market orientation, communication and transportation facili- 
ties, nature of product, local considerations like site 
value, attitude of local government including the planning 
restriction, labour availability, ancillary services, etc. 

1.1 Classification of Warehousing Problems ; 

Warehousing problems can be classified based on ware- 
house capacity, cost of opening the warehouse and customer’s 
demand characteristics. The variations of the problem in 
these three basic classes are given below; 

1 . Warehouse Capacity s 

Based on the capacity of warehouse, the problem is 
referred to as capacitated or uncapacited warehousing problem. 
In the capacitated problem, the warehouse capacity is assumed 
to be fixed while in uncapacited problem it is assumed to be 


infinite . 
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2. Cost of OTjening the Warehouse : 

If there is no fixed cost associated with the opening 
of the warehouse, then the problem is known as 'Depot location 
Problem'. On the other hand, there may be fixed costs of 
opening each of the warehouses. In such a case we have to 
determine the optimum number of warehouses and their loca- 
tions such that the overall cost of transportation and fixed 
charges of opening the warehouses is minimized. This is 
known as fixed charge warehouse problem. 

3. Consumer ' s Demand ; 

Depending upon the type of customer's demand the problem 
can be classified as Deterministic Demand Problem or Probabi- 
listic Demand Problem. 

The above stated classification scheme is depicted in 
Pigure 1. 

1,2 Motivation ? 

A considerable amount of work has been done by previous 
researchers on the fixed charge warehouse location problem. 
"Various analytical and heuristic approaches have been deve- 
loped for this problem when the demand is known. However, 
no attention seems to have been given to the case when 
demands are probabilistic. It was thus felt that structuring 
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the model to accommodate a pr ohahxlistic demand constraint 
would be worthwhile effort. Therefore, the present work is 
directed towards the optimal location of warehouses in a 
given system when demand is considered as probabilistic. 

The problem in its simplest form can be stated asi given a 
number of demand points where the demand from the customers 
are concentrated. The demands are assumed to be probabi- 
listic and are supplied from warehouses which are to be 
located at a given sets of sites. The objective is to 
determine the following such that the total cost which com- 
prises of transportation and fixed costs of operating the 
system is minimized; 

(a) The optimum number of warehouses to be opened. 

(b) The location for the warehouses to be opened (chosen 
from a given set of alternative sites) 

(c) The set of demand points to be served by a given 
warehouse. 

To seek logical answers to the above problem, an 
analytical model has been formulated. The stochastic nature 
of demand is accounted for by structuring the model as a 
chance constrained programming problem. The probabilistic 
constraints are reduced to deterministic constraints by 
using the normal distribu-tion. The basic philosophy of the 
branch and bound procedure has been utilized to solve the 



problem. The various branching rules have’been used to 
arrive at the optimal solution. 

A computer package for the solution methodology has 
been developed in f'ORTEJUM I¥ for the IBM 7044 system. The 
computer package is quite general in nature and in its 
present form can handle problems involving 50 warehouses 
and 100 customers. 

1.3 Organisation of the Thesis ? 

In addition to this introductory chapter, there are 
three more chapters. In Chapter II, a brief survey of the 
existing literature relevant to the kind of the problems 
with which this thesis is concerned j is given. 

Chapter III is divided into three parts. In part 1 
the problem is formulated. In part 2, solution methodology 
IS presented to solve the problem. In part 3 , an example 
IS considered to demonstrate the procedure. 

In the last chapter, computational experience with 
problems of different sizes is presented. Further extension 
of the problem and suggestions for future work are also 
described in this chapter. 



GHAPTBE II 


LITBiaTUIiE RBVIEV/ 

A review of the literature on location theory 
indicates that considerable research effort has been 
devoted to the fixed charge location problem. The problem 
in its simplest form involves the shipment of a single 
commodity to n customers with known demand from m potential 
warehouses. The warehouses which are opened# for making the 
shipments incur certain fixed' costs. The objective is to 
select from the potential warehouses, which should be 
opened to minimize the total cost, i.e. , transportation 
costs and fixed costs of operating system. 

The current state of art regarding the approaches 
available for solving this problem can broadly be cate- 
gorized as follows j 

(1) Heuristic Approaches, 

(2) Simulation and Graph Theoretic Tecnniques , 

(3) Branch and Bound Approaches, 

(4) network and Mathematical Methods. 

liisoussion of these approaches will be limited to 
some of the important contributions 
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2.1 Heuristic Approaches ; 

The difficulty encountered in working with heuristics 
IS that there is no way to know how far the answer obtained is 
from the optimal. This difficulty tends to decrease the 
validity of a sensitivity analysis by obscuring the cause 
of a change in the solution which could occur if a particular 
parameter is varied. Por example, results from a particular 
procedure may be same as optimal for cas® in which the 
fixed costs for operating a warehouse are large relative to 
the variable costs, but when that ratio is reduced the 
solutions produced move away from the optimal solutions. 

Since we may be interested in the sensitivity of the solutions 
to changes in this ratio as a part of the general sensitivity 
analysis, the conclusions from using the heuristic procedure 
will be confounded by the fact that the quality of solution 
IS changing as well as the solutions themselves. Inspite of the 
above difficulties, the heuristic procedures still offer 
several advantages. The advantage of speed is a major 
consideration and the comparisons with exact solutions have 
indicated that this is not gained at the expense of a large 
penalty in terms of quality of solutions. In addition, 
heuristic procedure are versatile, sometimes more so than 
, problem specific exact procedures such as we have in case 
of the warehouse location problem. This may mean that 
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additional factors or restrictions can te incorporated more 
easily into tiie lieuristic procedures than the exact proce- 
dures . 

Kuehn and Hamburger [12] have developed a heuristic 
or non-optimal algorithm to find good solutions to the 
warehousing problem. They assume that transport costs are 
linear with the amount shipped and facility costs are of the 
form, 

i'(y) = a+ by if the facility exists, 

1 1 11 

=0 if it does not exist. 

Thus F^(y 3 _) consists of fixed charge that is independent 
of storage and a linear cost which does depend on the amount 
of storage if the facility exists. Since the expansion cost 
b^ IS assumed for each facility, the demand points which a 
given facility will serve is determined simply by assigning 
a demand point to that facility for which the sum of shipment 
cost and expansion cost is minimum. 

The solution method begins with a single facility, 
and another facility is added to see if the total cost can 
be decreased. The assumption is that the best H facilities 
are contained in the set H + 1 facilities. Termination of 
solution occurs when it appears that another facility cahnot 
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be added without increasing the total cost, and those ware- 
houses which have become uneconomical are eliminated. A 
warehouse becomes uneconomical if a subsequently added facility 
can serve enough of its customers at lower cost. 

Peldman, Lehrer, and Ray [5] have assumed a more 
general form for the facility cost (continuous concave 
function). They assume, too, that the transport cost is 
linear in the amount shipped between points. Since 
does not have an expansion cost which is linear with the 
size of the facility, the problem of assigning demand points 
to existing facilities is much more difficult to accomplish 
than before. Whereas Kuehm and Hamburger begin with one 
facility, Feldman et al start with all facilities existing 
and drop out uneconomical warehouses. Termination occurs 
when no further savings can be achieved by elimination of 
any warehouse. 

Irysdale and Sandiford [3] have suggested a heuristic 
which occupies a middle ground between precise analytical 
modes and simulation. It uses a 'dropping heuristic' com- 
bined with steadily incrementing fixed warehouse cost as 
well as a local enumeration after each 'dropping iteration'. 
This approach offers certain advantages over rigorous methods, 
in terms of greater computational speed and flexibility in the 
form of the cost function. However, the procedure lacks the 



12 


ability to reach with certainty the precise optimum. 

2.2 Simulation and (rroup Theoretic techniques ; 

Simulation has been used by Markland [15] and Shycon 
and Maffei [19] for warehouse location problems . The pro- 
cedure essentially procedure a total cost curve as a function 
of the number of warehouses that are to be included in a 
physical distribution system. Samples of specific warehouses 
are chosen and to each warehouse those customers are assigned 
which can be supplied at the lowest variable cost and then 
total costs are evaluated. Ehumawala remarks [9] that in 
using the simulation approach, the need to solve the trans- 
portation programs increase the computation time for getting 
results . 

A group theoretic approach has been suggested by 
Kennington and Unger [10]. They have found that associated 
group problem can be especially structured with predictable 
inter-rela tionship among the group variables. The depen- 
dency of sub-group order on the route capacity is also 
described. Authors claim that by taking advantage of 
group theoretic properties of the problem the exact optimal 
solution can always be found. 
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2,3 Branch and Bound Appro8ch.es : 

The hranch and bound algorithm developed by little , 
et al [13] for travelling salesman problem was the starting 
point in the development of many branch and bound algorithms, 
for warehouse locational problem. Notable contribution have 
been made by Efroymson and Ray [7], Spielberg [18] and 
Khumawala [7>8]- ^1 these investigators have used the same 

philosophy for search in the use of branch and bound for 
solving warehouse locational problem. The efficiency of a 
branch and bound procedure depends mainly on the determina- 
tion of the lower bounds. The bounding procedures proposed 
are based on minimum total cost with integer restriction 
on the number of warehouses to be opened. 

The branch and bound procedure requires generating of 
a sequence of partial assignments and analyzing the comple- 
tion of each in search of successively better feasible solu- 
tions. Efroymson and ‘'Ray's method involves selective 
enumeration which is guided at each stage by a bound on the 
value of the objective function obtained at that stage. The 
efficiency of branch and bound m his formulation is due to 
the fact that non-integer result at each stage is an obvious 
solution to a simple linear programming problem. The 
algorithm terminates when there are no nodes (stages) left 
which can possibly produce an improved solution. A hic den 
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advantage in the formulation, in addition to the fact that 
it terminates optimally, is that the solution is independent 
of all linearities in the transportation cost function. 

Spielherg [18] has extended Efroymson and Ray’s branch 
and bound algorithm by incorporating features to attain 
better computational efficiency. In this paper, the author 
has reported his computational experience for varied sizes 
of problems. The largest problem considered involves 100 
warehouses and 150 customers. Spielberg has also suggested 
models to incorporate linear facility expansion cost and 
budget constraint on total expenditure. However, he has not 
reported any computational experience. 

Recently Khumawala [7] has developed various branch- 
ing rules. These rules are largest Omega rule, Smallest Omega 
rule, largest Delta rule, Smallest Delta rule, largest ’Y’ 
rule. Smallest 'Y’ rule , largest Demand rule and Smallest 
Demand rule. These rules reduce the size of the branch and 
bound decision tree and thus lead to considerable saving 
from computational standpoint. By using the information 
already available at each node through node simplification 
criteria the author has suggested an improvement in solving 
the linear program without any computational work. The 
computational experience and efficiencies of various branching 
rules have also been reported with large number of problems. 
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2,4 G-raphical and Mathcmetical Methods ; 

Marks [16] has formiilated a fixed charge prohlem 
which IS more general in nature. His model allows the 
facilities to be constrained in their capacity, and the 
warehouse is considered as an intermediate point between 
the sources of product and the demand points for the product. 
The facility cost function has a fixed charge and a linear 
expansion cost. If the source or demand points do not exist, 
the problem reduces to the general warehousing model. 

According to Mark's model, product is to be provided 
in specific quantities to each source area of finite supply. 
In mathematical form, the problem can be expressed as, 

m m n [E. £ 

Minimize; v F y + £ T. 

r=l ^ 1=1 d=l ^ i=l k=l 

Subject to the constraints; 
m 

^ k = 1,2,..., p 

i=l ^ ^ ^ 

h -, = ^ ^ =:l,r,..., m 

j=l ^ ^ k=l 

P 

r < Q^y^ 1 = 1,2,..., m 

kSL 

m , 

> y X« J = 1,2,..., n 

J — 1 J ~ 3 
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where. 


X_^ , Xt^ ^ are non-negative integer, y = (O, 1) 

1 J X ^ 1 


y^ = 1 if the i-th facility is built 

= 0 otherwise 


X^ = 
1 3 


k 1 


3 - 3 


R. 


flow of material from facility i to sink j 

flow of material from source k to 
intermediate point i 

0 + E = unit cost associated with a transfer 

3- J 

of material from facility i to sink j 

unit variable and ^associated with using 
sink j 




/ 

G, . 
ki 




1? 


S. 


= ^k "^1 ~ '^ 3 iit cost associated with 

transfer of material from source k to facility i 

= Unit shipping cost from source E to facility i 

= unit variable cost associated with using 
source k 

= Unit variable cost associated with using 
facility 1 

= Fixed charge for establishing facility i 

= amount supplied at source k 


= upper bound on amount demanded at sink 3 
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= lower bound on amount demand at sink j 
Q^_ = capacity of the i-th facility 

m = number of proposed facility sites 
n =: number of demand points 

p = number of supply points 

The solution technique is based on the recognition 
that a network may be applied to the problem. The graph re- 
presenting the network is shown in Pig. 2. A capacitated 
node for each facility has been added so that a capacity 
constfaint and a linear cost function may be ascribed to 
eacn facility. 

The method begins by approximating the fixed charge 
cost function with the linear unit cost 'Plie true 

cost function and the linear approximation cost fxmction 
are shown in Pig. 3. It should be noted that the approxi- 
mate cost function under -estimates the true cost function 
except when the flow through the facility is zero or Q^, 
the capacity of the facility. The approximation forms the 
basis for a branch and bound scheme. 

Assume all facilities are open and have the approxi- 
mate cost function. Solve this initial problem using an 
out-of-kilter algorithm. If the resulting solution is such 
that flow through each of the facilities is either zero or 
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the capacity of the facility, the optimal solution to the 
fixed charge problem has been found. If not, branching 
takes place, and a facility is either included in the 
solution set with its fixed charge added to the cost or 
excluded from the solution by setting its capacity to 
zero. This procedufo continues until an optimal solution 
IS found and verified. 

A recent and promising approach to the fixed charge 
problem has been presented by Walker and Lynn [20]. The 
thiust of their algorithm is to complement the evaluator of 
standard simplex procedure with the fixed charge associated 
with the entering vector and decrease the ovaluator by the 
fixed charge associated with the existing vector. In addition, 
the fixed charges associated with the basic variables driven 
positive in the exchange must be added to the evaluator, and 
the fixed charges associated with basic vectors driven to 
zero substracted. Let 

*P = vector currently not in the basis 
u 

Pk = vector currently in the basis which is to be 
removed if P enters, 

= fixed coat associated with any vector P^ 

X “a. 

6 = level at which vector P, is to enter chosen by 

J 0 

0 = Min lel (^' ) = X > 0, I the set of 

ij ko 


basis vectors 
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Zj-Cj = the conventional simplex evaluator 

Aj = the modified evaluator 

S = the set of basic variables vrhich are zero but 

become positive due to the exchanging and 

1 = the set of basic variables which were 

positive but become zero due to the exchange 

The modified evaluator is then 




- f. 


- 




i£b 


1ST 


In 8 minimizing problem, a negative indicates 

J 

an improvement in the objective function on the appropriate 

exchange. If iteration are repeated until all Ai are 

J 

positive, a local optimum will have been reached. Since 
better solutions may exist but not at adjacent extreme 
points the algorithm forces non-basis vectors into the basis 
even though their inclusion causes an increase in the 
objective. At the new extreme points, the algorithm again 
seeks favourable evaluators on the possibility that a new 
pattern of exchanges will lead to a better solution than the 
previous local optimum. 

To summarize, the literature survey reveals that 
various analytical and heuristic approaches are available 



for solvuig this prohlem whsn dSiiiard is known. However , 
no attention seems to have been given to the case when 
demands are probabilistic. Therefore present work is 
directed towards the optimal location of i^are houses in a 
given system when demand is considered as probabilistic. 
In Chapter III the formulation and solution procedure of 
this problem are given. 



GHAPTiiH III 


PROBLEM POPULATION AND SOLUTION PROCELUEE 


3 .1 Statement of the ProLlem i 

Ne are given a number of demand points where the 
demands from the customers are concentrated. The demands 
are assumed to be probabilistic and are supplied from 
warehousesa which are to be located at a given set of 
sites. The objective is to determine the following such 
that the total cost, which comprises of transportation 
and fixed costs of operating the system is minimized. 

(a) The optimum number of warehouses to be opened, 

(b) The location for the warehouse to be opened (chosen 
from a given set of alternative sites) 

(c) The set of demand points to be served by a given 
warehouse . 

3 . 2 Assumptions ; 

A mathematical model based on the following assump- 
tions has been developed. 



23 


(i) Unit transportation costs are constant, i.e,, they don't 
change with the volume transported and with time. 

(ii) Warehousing fixed costs are not affected by demand 
(ill) Each warehouse has unlimited capacity 

(iv) The demands are concentrated at a known set of points 

or locations. Demand points and warehouse location can coincide. 

(v) The demands are normally distributed random variables. 

The following notation is used for development of 
the mathematical model. 


3.3 Uotationi 


m 

n 




C 


ij 


^i 


number of possible warehouse sites 

total number of customers (demand points) 

fraction of the demand of customer i which 

IS satisfied by the warehouse located at 

site 1 = 1 ,...^ m, 3 = 1 ,..., n 

fixed charge resulting from establishing a 

warehouse at i 

unit transportation and handling cost from 
warehouse i to customer 3 

cost of supply the entire demand of customer 3 
from a warehouse located at site 1 , 

1, if a warehouse is located at site 1 
0, otherwise 
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P 


X 






n. 


r 


J 






a 





a(rj)= 


the set of customers that can he supplied by 
a warehouse at i , 

the set of warehouses that can supply 
customer j, 

the total number of demand points in the 
set 

demand of customer j (a normally distributed 
random variable), 

demand of customer j which is actually 
satisfied , 

minimum level of confidence with which 
customer's actual demand, i.e. r^ must be 
satisfied. 

abscissa associated with the left tail of the 

normal probabilitj’- density function for the 

confidence level , 

mean of the random variables r^ , 

standard deviation of the random variable r^ 


3 , 4 Problem Pormulation; 


Using the above notation the warehouse location 
problem stated can be formulated as follows: 
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Minimize 

Z = 

in 

f 

1=1 

n 

r c 

3=1 

X 

13 

m 

+ r E Y 

1=1 “ 

(1) 

S.t. 

r 

isH 

X 

= 1 

d 

— 1 ^2^09 *9 1^ 

(2) 


V 

X . 

< n.Y£ 

1 

2 ^ o 0 c ^ m 

(3) 


X 

> 0 


V 1 

and 3 

(4) 



= (0, 

1) 

V 1 


(5) 


P[D 





(6) 


Explanations for the logic used for the development; of the 
above stated mathematical model are given belox^; 

(i) The objective function gives the total cost of 
operating the system considering the cost of transportation 
and the fixed costs for the warehouses which are opened, 

(ii) Constraint (2) indicates that the demand for customer 
j IS satisfied by some combine. tion of supplies from the 
warehouses , 

(ill) Constraint (5) indicates that the sum of inactions 
of satisfied demands of the various customers supplied by 
the open warehouse located at site i cannot exceed the total 
number of customers. However, this sum will be zero when 


the warehouse is closed . 
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Oonstraint (4) represents a non-negativx ty restric- 
tion on the decision variables X . X will be zero when 

^ J 1 J 

warehouse is closed. 

(v) Bq^uation (5) represents e zero-one type of res- 
triction. IS 1, if the i-th warehouse is open, otherwise 

IS zero, 

(vi) Constraint (6) gives the level of confidence with 
which the j-bh customer's demand will be satisfied. 

It IS interesting to point out that if there were 
no fixed costs associated with the opening of warehouses the 
optimum solution to the problem would be locate a warehouse 
at every site. In such a case the problem reduces to a 
simple assignment problem and one is required to determine 
the warehouses which should supply a particular customer. 

On the other hand, if the transportation costs are ignored 
the best solution to the problem would be to locate one 
warehouse at the site with the smallest fixed cost. Thus 
we see that in the general problem where both the tra.nspor- 
tation costs and the fixed costs are considered a balance 
between these two types of costs need to sought so as to 
minimize the total cost of operating the system. 

The mathematical formulation as given by (l) to (6) 
carries the stiucture of a chance-constrained programming 
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problem because the demand constraints given by (6) are 
probabilistic. The chance -const rained programming problem 
can be converted into deterministic problem using the 
following approach. 

5.4.1 Reduction of Cha ncc-g cnstrcine d Prograrm ein c Treble : 
to A Deterministic Problem; - - - 


let r^ be a normally distributed random variable with 
|i(r ) and o the mean and standard deviation of r , respec- 
tively. The parameters u,(rj, cr are assumed to be 
estimated from past experience. Then, 

P 


P [1^ - 

orT ^ 


|a(r ) 


or. 


] > a. 


The terms on the left hand side of each of the brackets- d 
inequalities are standardized random variables with zero 
mean and unitary standard deviations. 

let Z be the abscissa associated with the left 

a j 

tail of the standardized normal distribution function for 
probability level a , as illustrated in Pig. 4. 

J 
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Then, 


r[S aj ^ -i 




ar. 


] = 


a 


ar 


> X or 
- aj 


T) > X„, a, + Li(r.) 

3 - aj r3 ^^3^ 


Therefore for a given choice of the decision variable, 1 , 

J 

we can write 

xff + ,(r^) (7) 

Since there is no advantage in shipping more than the 
required amount, constraint (7) can be written as, 


Substituting relation (7) for relation (6), the warehouse 
location and shipping problem can be written in its deter- 
ministic equivalent form without recourse to a probability 
for the constraint. The mathematical formulation of the stated 
problem reduces to, 
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in n 

Minimize ^ = Z H err + |j,(r )) t X, 

T tHi J j J 1 3 ■ 


m 


+ z 
1=1 

P 

1 

Y 

1 


(9) 

z 

ISX 

J 

X 

13 

= 1 J = 

1,2,.. . ,n 

(10) 

£ 

X 

nY 1 = 

1,2,. . . ,m 

(11) 


1 J 

~ 1 


^ n > 

IJ - 

0 

i and 

J 

(12) 

^i = 

(0, 

, 1) ^ i 


(13) 


.4.2 Wecessery Conditions of Optimality ; 

From Appendix I (equations (18) to (22)) and equation 
7), we have following relations, 

A ^ n g. 

Minimize Z = n T 

ieK^UK2 a=l ^1 

(14) 

X = 1 J = 1, 2,. . . , n (15) 

J 

t ^ 1 (16) 

3 =1,2,..., n 


( 17 ) 
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The legrangian espression for the problem can be given as 


n 


\j = V r (I, ^ t , + -^) i 

J 5 ieKj_TO 2 J=1 ^ ij n ij 


n 


' L - ""b °b - 

J 


>. 5 , 


V tij X - > >0 

i^^UK2 ^ 


At the opto-num value when we have decided bo open the i-^s 
warehouse we have, 


h*: - ^ ° 


(18) 


If A IS interpreted as the imputed value of a delivered 

t] 

item, then (14) indicates that in the optimum case the 

cost of shipping one unit to destina.tion j from warehouse 

1 , t must be greater than or equal to the imputed 

^ 3 

Value of a delivered unit, 

It should be noted that Xj represent the change in 
the minimum attainable cost resulting with tne reduction or 
one unit in the quantity [Xa^ ar^ + whose value 



depends upon the selected values of a and and the 

estimates of the parameters or and |i.r , of the probability 

J J 

distribution. 

Suppose that the manager of a shipping firm has 

signed a contract with the customer which demands that the 

shipping firm must meet the actual requirements of the 

customer 50 percent of the number times the order for 

shipment is placed. In rhis case a «= .5, and = 0. 

j 

The constraint (15) would involve now only the expected 
values of the requirements distribution, pCr^). In short, 
the delivery firm is only concerned with the first 
moment of the requirement distribution. In this case > 

J 

measures the change in tne minimum attainable cost 
associated with a 'small' change in !i(r ). 

J 

If the shipping firm manager feels that the 3-th 
customer's actual requirement must be met more than fifty 
percent of time, then a will be greater higher than 
0.50 which implies thab Xq- will be greater than zero. 

J 

For given or ) end a , the effective requirement 

J J J 

for customer 3 is [Xq., or + pCr )]. The quantity 

J J J 

Xq. or can be thought of as an additional quantity that 

J fJ 

must b© shipped (over expected requirement i.e. when 
a = .5). A decrease in a can be determined so that 

J 
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tiie value of crr^ is reduced ly one unit. The effect 

of such a change in a on the ninimum attainable cost is 

x) 

measured by A . let Xq. correspond to the value of 
J J 

which results in a unit reduction in the value of 
[Xa crr^ + ii(r )]. Then o' - a- = where a' 

J O J J J J J 

corresponds to X' 

3.4.3 SelectiOxi of a : 

J 

The value of a to be used in expression ( 8 ) is 

J 

normally treated as s management decision parameter. The 
management can decide the value of cc either using pui-oly 

J 

intuitive approaches or by coupling the intuitive approaches 
with the systematic procedure given in the following 
paragraphs. 

The selection of the minimum confidence level should 
be based on a balancing of two different ’costs*. On one 
hand tiers is the cost of violating a constraint, i.e., of 
not meeting the minimum requirement at bhe destination. 

These 'costs' must be determined outside the model. In 
general, the higher the level of a selected the lower the 

J 

expected level of these costs. On the other hand there 
are costs ineasured in terms of the objective function of the 
model which increase as the confidence levels increase 
since the problem then becomes more severely constrained. 
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The latter cost can be estimated using parametric programming. 

It would be worth while if both costs could be captured in the 

objective function by including a specific cost associated 

with the confidence level with whicb a constraint will be 

violated. In short, tne determination of levels for a is 

properly apart of the optimizing model. This technique 

require a means of weighting the costs given in ( 9 ) and other 

costs associated with the selection of a . The weights 

d 

must be determined through the use of a utility function 
for the firm's management. 

If we let u (a ) be the 'cost' associated with the 

J 

violation of the j-th constraint then firm' s utility 
function can be written as UfZ, 'j( q:2^ • • • 5 

If we fix each of the n's and minimize Z we have found 

cJ 

one point on the efficiency frontier for this firm. The 
complete frontier can be obtained bj?- paraure trie ally altering 
the a 's to find further best points on the frontier, 
of the firm. 

The costs associated with not meeting a customer's 
minimum requirements will undoubtedly be difficult for the 
firm to measure since the loss incurred majr not bo purely 
monetatry. Instead, the firm is more likely to be concerned with 
some estimates of frequency wit!a which a customer's require- 
ment will not be met even though it may not be able to 
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convert the freouency directly into cost. In other vords, 
the manager of a deliverv firm maj'' not he ahle to estimate 
the cost of not satisfying the requirement of a particular 
customer. However, he may have some subjective feelings 
about how often, on the average, he can disappoint each 
customer without dire consequences. In snort, the a^. 
parameter used in relation (6) may be taken to be the sub- 
jective estimate of this critical quality of service. 

Besides, the utility function approach, the value of 

a can be determined by considering the value of number of 
d 

units lost. Let UI represent the number of units lost for 

the j-th customer when the shipment is planned accordixig to 

the desired a . Mathematic aHly UL is given by, 

J J 

CO 

\ (r -d ) f(j6) dr 
d J J ^ 

J 


Since, the demands are assumed to be normallj?' distributed, 


UL 




1 


J . 'U - L’ «p - 2 ^ ' -3 


r - ir(r ) 2 
i (...i - ^) dr. 


rhr 


. ^(r -u(r ) 2 

3 

rrn -j "j -1*“!'^ ( ^ -I ) 2. 


= dr 


1 


J 


exp 


, d -ii(r.,) 2 

1 + (u(r^)~d,Ml- 


■oc 


C_J 
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The management cen convert tne number of units lost 

into monetarjr value, ML-,> by assigning appropriate weight 

J 

per unit lost for each customer, 
i'lL = V-f X UL 

V J J 

The value of ML decreases as a increases. On the 

J J 

other hand, the value of tiie objective function increases 
with increase in the value of a . Therefore, the optimum 

J 

corresponds to the point where the values of these two 
functions are equal. 


3 «5 Solution Methodology ; 

The problem as formulated in (9) to (13) is first 
solved as a linear-programming problem, igroring the integer 
restriction on Y^. If all the Y^ values obtained b^ solving 
the linear program are integers, then the problem has been 
solved. The solution procedure for solving the linear program 
IS given below. 


The assumption of unlimited warehouse capacity ^ives 
the linear program a structure which can be exploited to 
advantage as shown in Appendix I. Referring to Appendix I, 
we observe that, 



1 If C 


ij 


s 


n 


Min [0, + — ] 


0 otherwise 


( 20 ) 
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and 

= 0 , for 1 £ Kq 
= f Y. /n , for 1 £ 

feP^ ^ ^ ^ 

= 1, for 1 £ 

where , 

~ ^k' ^ ^2 

= 0, for K £ 

If for some K, is non-integer, then problem can 

be solved using number of methods available in the literature, 
imong these methods, the important ones ares enumerative tech- 
nuques by Efroymson and Ray [4]. Khumawala [7,8] and 
Spielberg [18] , group theoretic techniques by Kennington and 
Unger[10]; and heuristic procedures by Reldman, Lehrer and 
Ray [5j end Kuehn and Hamburger [12] . 

The problem of determination of the optimum number 

of warehouses selected for opening, out of the given m 

potential warehouses, is highly combinatorial in nature 

and the number of combinations which need to be considered 

increases as m increases. The total number of combinations 

m 

involved would be 2 -*■1. I>ue to combinatorial nature of the 
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problem, mostly computer oriented algorithms are developed 
to seek the optimal solution(s). The efficiency of a 
solution algorithm is generally characterized by the compu- 
tational time and storage requirements. Khumawala [7] 
end Efrymson and Ray [4j claim that branch and bound 
procedures can be utilized efficiently for this class of 
problems. In the following sections the basic philosophy 
of branch and bound procedure is applied to the problem 
under investigation. 

3.5*1 The Branch and Bound Algorithm ; 

The branch and bound procedure requires generating 
a sequence of partial assignments and analyzing the comple- 
tions of each in search of successive better feasible 
solutions. There are many systematic ways of choosing a 
sequence of partial assignments that provides optimal or 
near optimal solution(s). The procedufe we have employed 
can be represented by the flow diagram given in Rig. 5 . 

At the start of the algorithm, all warehouses are 
initialized as free and various node simplification 
criteria are applied. If at any node a free warehouse 
exists then branching rules are applied for fixed closing 
or opening of that warehouse. This procedure of branching 
and bounding continues till a node is reached where all 





K.f fc-H-Ki* i 

^hr, ;ca'?’'ii':r'..;ic': 


SOIjVK l 




t'Daik? : 


.>tBi‘£XC 
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values are integer. After re pching the node where all Y^'^s 
are integer, a linear program is solved. The procedure is 
repeated with all branching rules and comparison done 
among the solutions obtaiiied by various branching rules. 

The rule(s) which give(s) least objective function value 
IS the optimal solution. 

Further discussion on branch and bound algorithm 
requires familiarity with the definitions given in the nejct 
section. 

3. 5. 1.1 A Few Definitions ; 

Terminal Node ; A node where all values are integer is 

called a terminal node. 

Non-Q^erminal Node: A node which involve at least one 

fractional Y^ value is known as non- 
terminal node. 

Dangling Node: A node which may produce an improved solution 

' IS callwd a dangling node. 

Closed Branch: The branch along which the warehouse is 

constrhined closed is called closed 
branch. 

Open Brandh: The branch along which the warehouse is cons- 

trained open IS called open branch. 
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3.5.1. 2 bimplifxcations at Node: The following three node 
simplification criteria developed "by Jf roj'Txison and Ray [4] 
have been used for reducing the number of branches at each 
node. These criteria are applied sequentially. 

Criterion I; Ror Rixed Opening of b^arehouse : 

This node simplification criterion determines a 
minimum bound for cost reduction in 'opening’ a warehouse. 
If this bound is positive, the warehouse is fixed 'open’. 

Mathematically , 


For 1 e and 3 e 


V 


13 


Min [ Max ( 0^ • 

KeNj O (S^UICg) 

K^i 




0)J 


( 21 ) 




leP, 


- F 

ij 1 


( 22 ) 


Ifhi^ > 0, then = 1 for all branches emanating 


from that node. 


-v measures the minimum cost saving for customer 3 
X 3 

that can be made if warehouse 1 is opened considering all the 
non-closed warehouses at that node. Clearly if the sum of 
such minimum savings for warehouse 1 over all customers that 
it can supply exceeds the fixed costs, it will always 

pay to open warehouse i at this node. 
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Criterion II. Por Reduction of n^: 

This criterion reduces the total number of customers 
■which Can he supplied by a particular warehouse. 

Mathemati cally , 

For 1 e K 2 and j e 

If 

Min ~ ^ 10 ^ - ° 

(23) 

then is reduced by one. Of course, if the inequa,lit 7 
holds for all j e P^, the P^ = 0. Therefore P^ = 0 anc 
= 0 for all branches emanating from that node. 

Clearly, if an already open warehouse can supply s 
customer j cheaper (in terms of lower variable costs) then 
any of the ’free' warehouses at that node, then such a customer 
should be supplied by the open warehouse. Therefore, such a 
customer should not be considered as a potential customer 
for the free warehouse at that node. 

Criterion III; For fixed Glosina, of '’•varehouse ; 

This simplification deteimines a maximum bound on 
the cost reduction for opening a warehouse. If this bound 
is negative, the warehouse will be fixed closed . For 
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X e ^2 s-^ad 3 e the net saving, , for the customer 3 

when supplied from warehouse 1 can he expressed as, 

oO,. = Min [Max (G, -C_, 0 )] (24) 

KeM HK 

J 

The total savings 5 -.'2 which accrue from opening the 
warehouse 1 will be, 




(25) 


If — '’^3_ ^ 0? then = 0 for all branches emanating from the 
node under consideration. It needs to be pointed out that'v^,^ 
is similar to except th8.t for comparisons 

are made only overall the fixed open warehouses and u :■ is 

1 J 

the minimum saving for supplying customer 3 from the 
opened warehouse i. Clearly, if the sum of such savings for 
warehouse 1 overall customers that it can supply fails to 
exceed its fixed cost P^, then such a warehouse should be 
closed and eliminated from further consideration. 


for each node these three node simplification 
criteria are applied sequentially till no lurther simplifi- 
cation is possible. It should be noted that if criterion I 
fails then Criterion II and III cannot be applied. The cycle 
through simplification is an important part of the procedure 
hence it is described by means of a separate flow chart 


A'T veri T n T?! a- A 




Ql£j 


-'ij' o rwM 










45 


5 . 5 - 1.3 Branching Euless 

After applying all the simplif icabion rules at 8 node 
if K2 7!^ 0 then a warehouse is selected from the set of free 
warehouse at that node for further branching. The selected 
warehouse is either constrained closed or constrained open 
to yield an additional node. Instead of selecting the 
warehouse randomly j the following decision rules developed 
by -Thumawala [8j are used to arrive at an appropriate 
node for consideration. 

De lta Eule ss 

In the node simplification criterion I, 1 ^'s are 

computed for each free warehouse at every node. If h; ^ > 0 

then warehouse 1 is fixed open for all branches emanating 
♦ 

from that node. However, for the warehouses for which s 
are negative no decision is taken. Warehouses having 
large (i.e. small negative)'^ values are likely to be 
opened in the terminal solution from this node. On the 
other hand 5 those warehouses having small (i.e. high negative) 
devalues are likely to be closed in the terminal solution 
from this node. We have two branching rules based on ’s. 

(1) j-jargest Delta Rules Select the free warehouse for fixed 
opening which has the largest from the set of free 
warehouses. 
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(2) Smellesi; Delta Rule; Select the free warehouse for 
fixed closing which has the smallest f roi the set of free 

warehouses. 

Omega Rules ; 

In the node simplification criterion IIIj— h^'s are 
computed for each warehouses at each node. If _i^ <, 0 then 
warehouse i is fixed closed for all branches sma,nating 
from that node. however, for the warehouses whose h- are 
positive no decision is taken. Therefore vrarehouses having 
largest positive -h- values are likely to be open in the 
terminal solution reached from this node and vice versa. 

Two branching rules based on considerations are ; 

(1) largest Omega Rule; Select the free warehouse for 
fixed opening which has la^rgest positive n from the set of 
free warehouses. 

(2) Smallest Omega Rule; Select the free warehouse for 
fixed closing which has the smallest positive k'. from the set 
of free warehouses. 

Y Rule s ; For applying these rules a linear programming 
problem is solved for the selected node. Whether a warehouse 
IS fixed open, fixed closed or free depends on the value of 
Y^. If Y^ = the warehouse is fixed open, if Y^ = 0, the 
warehouse is fixed closed and if Y^ is fractional, the 
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■warehouse is free. A free warehouse whose Y is close to 
one is more likelv to he open in the terminal solution 
reached from that node as compered to the warehouses whose 
Y IS less. Con-versely the warehouse -whese Y is close to 
zero IS likely to he closed in the terminal solution 
reached from the node. Tnis leads to two branching rules 
based on the Y' s . 

(1) Largest Y Rule; Select the free warehouse for fixed 
opening with the largest Y from the set of free warehouses 
at the node having fractional Y. 

(2) Smallest Y Rule; Select the free warehouse for fixed 
closing with the smallest Y from the set of free warehouses 
at the node having fractional Y. 

Remand Rule ; The rationale here is that if a warehouse is 
supplying 8. large demand currently it is likely to be open 
in the terminal solution reached from the node and vice 
versa. The two branching rules based on demand considerations 
are ; 

(1) Largest Demand Rule; Select the free warehouse for 
fixed opening which is supplying the largest demand from 
free warehouses at that node. 

(2) Smallest Demand Riile; Select the free warehouse for 
fixed closing which is supplying the smallest demand from 
free warehouses at that node. 
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3. 5 •1.4 Use of brenciiing Rules for Node Selection. 

TRe use of the brench and bound procedure involves 
large computer storage and high computational time. This is 
because the number of nodes that need to be evaluated to 
reach an optimal solution has an exponential relationship 
with the total number of integer variables in the problem. 

For reducing the number of nodes which need to be considered 
for evaluation, node selection procedures developed by 
Khumawala [8] have been used. The node selection procedure 
IS used in conounction with each branching r\iLe. In essence, 
the node selection procedure can be stated as, 'Trace the 
preferred path which would lesult from an application of the 
particular branching decision rule'. 

The word 'preference' is used here to mean the 
motivation on which the free warehouse is selected. As the 
branching nle is applied in the branch and bound algorithm 
from the initial node [node 1, Fig. 7] we ex;pect the best 
solution to be along some preferred branch. If the node 
resulting from the preferred branch is not terminal 
branching will continue along the preferred path till 


a terminal node is reached. In figure 7, the preferred 
path will be l-3-5-7“ n. 
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Illustrative tixample I; (I'or Selection of a ) 

3 

A numerical example is considered here to illustrate 
the methodology for selection of cc.. 

tJ 

Conforming to the notation, the input data is: 

m = 10, n = 15, p(r ) = 26.00 for 'tf and - 5 t3 

J do 

and or = 5.3 for all V 
J J 

The fixed costs for the various warehouses (1,2, ...10) 
are (60, 185, 95, 65, 120, 70, 150, 90, 50, 80) respectively. 

The values of (■tj_j) are given in Table 3.1. It is 

assumed that for all columns a = a a given constant. The 

J 

various a's considered are .5, .6, .7, .8, .9 and .95. 

For various values of a, the values of satisfied 

demand d , the value of units lost (UL ), the value of the 
J J 

objective function are presented in Table 3.2. Fig. 8 is a 
graphic representation of total cost and number of units lost 
for various values of a. By assigning appropriate weightages 
to the unit lost, the unit lost curve can be converted into 
a cost curve. The super imposition of these two cost curves 
will result in a curve with a unique minima. The optimal 
value of a will correspond to the point where the super- 
imposed cost curve has a minimum value. 



Table 3.lJ Unit Cost Prom Warehouse i to Customer 
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Ta^le 3.2 


a 



Z 

.5 

26.00 

31.80 

18900 

. 6 

27.3430 

21.34050 

19892 

.7 

28.7804 

14.93820 

20921 

.8 

30.4621 

8.96370 

22125 

.9 

32.9751 

3.21735 

23796 

. 95 

34.7227 

3.042975 

25176 



X 10^ VALUE OF OBJECTIVE FUNCTION 



NUMBER OF UNITS 
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Illustrative Example II. (Eor the Procedtrre) 

A problem of (4 x 6) has been structured and solved 
to illustrate the various steps. The input data for the 
problem is given in Tables 3.3 and 3.4. 

Table 3.3 


Customers 1 2 3 4 5 6 


Desired confidence 
level for meeting 
the demand 

Mean Demand 
M.(rp 

Standard 

Deviation 

J 


Table 3.4; Unit Cost Matrix (4^^) and fixed Cost 


Ware- 

house 

Fixed 

Cost 

1 

Customs rs 

2 3 

4 

5 

6 

1 

84 

1. 2 

9.0 

7.0 

9.0 

14.0 

4.5 

2 

60 

5.0 

8.0 

6.0 

16.0 

12.0 

12.0 

5 

120 

1.1 

10.0 

8.0 

10.0 

15.0 

15.0 

4 

72 

1.3 

12.0 

9.0 

7.0 

11.0 

8.0 

step I; 

Find 


1 . e , the 

absicissa 

associated with left 


tail of the standard normal distribution function for the 
confidence level cc . Table 3.5 represents the effective 

J 


.8 .87 .9 .85 .93 .95 

12.0 5.0 16.0 10.0 25.0 30.0 

.4 .01 .1 .54 .8 2.1 
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demand of each, customer which is obtained by the expressiono 

Table 3.5 


Customers 

'‘a, 

1 2 

.8418 .1126 

12.34 5.01 

3 

.1282 

16.12 

4 5 

.1036 .2055 . 

10.55 26.64 3 

6 

1645 

3.45 

step 

II 

Obtain (O^j) matrix by multiplying demand 

of 

J-th 

customer with corresponding 

column of t matrix, 

-L J 



Table 3.6: 0 Matrix 





Customer 



Warehouse 

1 2 

3 

4 5 

6 

1 


14 .80 45.10 

112,89 

95.03 373,02 

150.55 

2 


37.01 40.09 

96,76 

163. 95 319,73 

410,74 

3 


13.57 50.11 

129,02 

105,59 399.66 

501.34 

4 


16.03 60.13 

145.15 

57.90 293.09 

257. e4 

Step 

III; 

Solve the problem as linear pro gram mi /ig 

problem 


Without integer restriction on Y^. from relation (20) the 
value of +(g^/n^) ; are calculated and given in Table 3.7. 

If all the are integer tnen problem is solved otherwise 

go to Step XT', 






Table 3.7; 

G + 

ix 

n 

1 



Ware- 



Customers 



house 

1 

2 

3 

4 


6 

1 

28,80 

59.10 

127.89 

109.03 

387.02 

164.55 

Q 

47.01 

50.09 

106 . 76 

178.95 

32-J.13 

420.74 

7? 

J 

33.50 

70.11 

149.02 

125 . 59 

419.66 

521.84 

4 

28.03 

72.13 

157.15 

85.09 

305.09 

279,64 


Therefore X^^ 

= X 22 = 

: ^23 = X 41 


= ^/ c =1 

45 

and all 

other 

ho =« 

and Y^^ 

= 1/6, Y 2 

= £/6, 

Y 3 = 0 /b 


and = 3/6. 

Since all Y^ are not integer therefore, procedure will 
continue. 

Step IV; Initialize the node and all war-houses as free. 

Step V; If it is the first node go to Step VI, if any 

other and the hranch is open go to Step VIII otherwise to 
Step XII. 

Step VI; Galculate and Check if any ^ ^ is 

positive, open that warehouse otherwise go to Step XIV, 


■^11 “ *^’^2 ~ ^13 ~ 

^21 ~ 6> ,^22 ~5«01j ^^=lo.l3' 24 ’ '' 
=1.23,^32=0^ "^33=0^ '^ 34 - = 0,1/33 
V 41 = 0 , \7^ 2=0, ’1'43 = 0 , "^4^ -21 . 07 


■^5 = 0 , -2^ = 0 

= 0 - ^-'^36 = 0 
, V^5=26.64,^3=0. 
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How 5 


^'1 

= 117.09 ~ 8n 

= 33.09 

'^2 

= 16.13 + 5.01 

- 50 = -38.86 


= 1.23 - 120 = 

-113.76 

<1 

= 21.07 + 26.61 

- 72 = -24. 29 


Because 0, hence vTarehousc 1 is fixed open, 

therefore = [l] and =- [2,3,4]. 

Step VII j If ^2 = go to Step XVII. 

Step VIII: If 8 particular customer is best supplied by ac 
already fixed open warehouse it is removed from further 
consideration. Here customer 6 is being supplied the cheapest 

Table 3.8 


Warehouse 

1 

Customers 

2 3 

4 

5 

6 

1 

14.80 

45.10 

112.89 

95.03 

373.02 

150.55 

2 

37.01 

40.09 

96.76 

163. 95 

319.73 


f-7* 

D 

13.57 

50.11 

129.02 

105. 59 

399.76 

- 

4 

16.03 

60.13 

145.15 

73.90 

293,09 

- 


Step IZ: Update K 29 if 1^2 ~ s S® Step ZVIT. 

Stop X : Calculate lu ^ and-C. for free warehouses. If any 

ij 1 

is negative, close that warehouse otherwise go to Step XIV 


1 
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5.01 + 16.15 + 55.29 - 60 = 1A.’3 

1. 25 ” 120 = -118. 77 
21.15 + 79.95 - 72 = 19.06 


Here-'- j 0 hencG it will befixed close and thus 

Kq = [3], = [1], Zg = [2, 4]. 

Step XI i If go to Step XTII, 

Step XII; Delete the customer which, can be best supplied from 
an open warehouse forthc further consideration for free 
warehouse. liere customer 1 is being the best supplied after 
closing of warehouses 5. 


Table 5*9 


Warehouse 

1 

Customer 

2 5 

4 

5 

6 

1 

o 

00 

H 

45.10 

112.89 

95.05 

575.00 

150.55 

2 

•T 

- 

40.09 

96.76 

168.95 

519.75 

- 

; 

4 

- 

60.15 

145.15 

73.90 

295.09 
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btop ZIII; G-o to Step IV. 

Stop ..I Vs If K 2 = [0] goto s-cep I.VII. 

Step XI s Select the warehouse according to tli^^ branching 
rule that is being applied, 1 . c. if largest delta rule is 
being applied then warehouse 4 will be fixed open since 
IS the largest among all h^.'s. 

hs 2 ~ "58.86 

= -24.25 

Step XVI; G-o to Step V. 


Step XVII; Solve linear programming and wc find 

X-, 1 = 1 , L. ^ = 1 , X, ^ = 1 5 X ^ = 1 , X * .j — It 


all other X = 0 and value of objective function is 

1 J 

846 . 565 . 


Step XVIII; Go to Step IV until all branching rulc° are 
exhausted . 


Step XIX; Compare various solutions obtained by the difforoiit 
branching rules and obtain minimum value. 



GHAPTLii IV 


BESJLTb AjTD PlbOUSSIONS 

A computer prograic of the methodology’’ proposed in 

Chapter III has been de-veloped. The program is written na 

fortran IV. A listing of the program is given in Appendix III. 

The input data requires the values of m, n, r^j u(r^); ^(r^), 

t , F . For the validation of the procedure a problem from 
1 J 1 

literature given by EIhumax:ala [8] has been selected. Since 
this problem assumes deterministic customers' demand, a 
value of a = .5 has been used for all the customers. The 
results obtained using the proposed methodology are found 
to be the same as those of Ehumawala. For evalua.ting the 
computational efficiency and limitations of the solution 
methodology, twenty problems were generated and solved. In 
all, problems of 4 different sizes, viz. (10x15), (153^20), 
(20x25) and (25x50) have been considered. The first number 
corresponds to the number of warenouses and the second 
number represents the number of customers. The computational 
time requirement for a problem of given size has been 
determined based on the average time required for solving a 
set of 5 problems of the same size. In Table 4.1, computational 
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tiine reauirements for problems of venous 

sizes 

are tabulated. 

Ixi Table 

4. 2, 

the a 

iverage 

times for using 

rhe various branch- 

mg rales 

IS given. 

The 

perfomance of various 

branching 

rules to 

generate optimal 

solutions IS tabulated 

in Table ^.5. 

Table 4 

.12 

Average Computation Time Required Por 
Problems of Different Sizes. 

ilo. of 

Vare house 

No. 

of Customers 

Average ^me in secs. 


10 



15 


20 


15 



20 


41 


20 



25 


44 


25 



50 


55 

Table 4 

. 2: 

Computational Time Requirements for 

Various Branching Rules. 

Size of 
Problem 



Branching Rules 




LDR 

LOR 

LYR 

IDR SDR SOR 

SfR 

SDR 

10 X 15 

. 2 

.4 

. 6 

.2 .2 .4 

.6 

.2 

15 X 20 

1.0 

1. 2 

1.5 

1.0 1.0 1.5 

1 . 0 

1.0 

20 X 25 

1.5 

1.7 

2.0 

1.9 1.7 1.0 

1.8 

1.0 

25 X 50 

6 

7 

12 

6 6 8 

20 

6 
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TeTole 4.3 5 Performa/ice iveluation of Various 
Branching Buies. 

(frequency of Obtaining Optimal Solution 
with Various Branching Buies out of 5) 


Problem 

size 

LDR 

LOR 

LYR 

LDR 

SDR 

SOR 

SYR 

SDR 

10 X 15 

5 

5 

4 

5 

5 

5 

5 

5 

15 X 20 

5 

5 

4 

5 

5 

5 

5 

4 

20 X 25 

5 

5 

4 

4 

4 

5 

4 

4 

25 X 50 

5 

5 

2 

3 

4 

3 


2 


The following important inferences have been drawn 
from the results. 

(1) Regardless of the branching rule applied, the 
computer times required for solving the 4 different sizes 
of problems considered, is very low (Refer Table 4.1). 
Therefore, the procedure can be used for fairly large 
realistic warehouse location problems with reasonable 
computer time. 

(2) It IS observed that for problems of smaller size, e*g., 
(10 X 15, 15 X 20), all branching ruies give the optimal 
solution but for large size of problems e.g. (20 x 25, 

25 X 50), Largest Delta and Largest Omega rules generaU.y 
give the optimal solution. 
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(3) I rules always take .-tiore time. This can be attributed 
to the rfact that applicaticn of these rules need solving of 
linear programming. 

(4) Jjince Largest Lelta and Largest Omega rules generally 
give optimal solution, it is, therefore, recommended that 
whenever one is reluctant to apply all these rules to arrive 
at the optimal solution, these two rules should be applrsd. 

I rules should be given least priority since they always 
take more computational time. 

bug^estions for future Work ; 

In the beginning, whenever a real life situation is 
modelled the model has to be, of necessity, kept simple 
either because the actual constraints cannot be included in 
a simple way or because solving the model becomes intractable. 
To mirror real life processes the model can be improved and 
extended in ma^ny wa3''S. Below wc mention a few directions in 
which further work can be done towards this end. 

1, The model presented in this thesis i& based on the 

assumption that the customer’s demand is nozmally distributed. 
Instead, the demand could follow other general distributions 
like Poisson, Binomial and Exponential distributions, j.t 
would be interesting to develop models and the solution 
methodologies for situations when all the customer’s demand 
don't follow the same distribution. 
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2. In the present model the a 's are being determined 

J 

outside the model. However, a model needs to be developed 
where the de toriiination of optimal a is treated an integral 

J 

part of the model. 

3. It IS assumed tha 'c warehouses are uncapacitared. 
However, in real life situation, there is always a limit 
on the size of the warehouse. Therefore, models for the 
location of capacitated warehouses considering probabilistic 
demands need to be developed. 

4. Hor an exercise like tins to be more meaningful, 

the system must be considered an a whole. Goods are rarely 
transported directlj?' from plant sites to customers, that is 
to say that warehouses are not located at the plant. Normally 
goods are shipped from the factory to the warehouses and 
thence onward to the customers. This two tier transportation 
model should be included in the formulation for a more 
realistic approach. 

5. The model neeefe to be furtcer extended for nonlinear 
transportation costs and constrainis like budget, etc, 

6. In implicit assumption in this work is that we are 
working in a discrete space. However, it would be worth- 
while to develop models for locating warehouses in 


continuous space. 
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7. In tins tiiesis, eight heuristic branching rules have 
been used and evaluated, however, one could conceive oH 
many more branching rules. A few more rules which can be 
tested are: 

Largest n^ Rule: At a given node, from the set of free 
warehouses, select the warehouse which is supplying largest 
number of customers for fixed opening. 

bmallest n^ rule; At a given xiode, from the set of free 
Warehouses, select the warehouse which is supplying smaillest 
number of customers for fixed closing. 

Largest Average Cost Rules At a given node, from the set 

of free warehouses , select the warehouse which is supplying 

n 

at the largest average cost ( ^ (C +? )/n ) for fixed 

J =1 ^ 

closing. 

bmallest Average Cost Rules At a given node, from the set 

of free warehouse, select the warehouse which is supplying 

n 

at the smallest average cost ( ^(C +f )/n ) for fixed 

J=1 ^ ^ 

opening. 

8. The node simplification procedure used is based on 
simple logic. xTowever, some mathematical approecnes for 
node simplification need to be developed. Hopefully, a more 
sophisticated node simpl t? '-eetior prooeduro will reduce the 


computational effoiC- 
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APPENDIX I 


SIMPLIFIGAIION OP LIKBiR PROGRAIi 

The assumption of unlimited warehouse capacity gives 
the linear program a structure which can he exploited 
advantage as shown helow: 

¥e have, 

m n m 

Minimize Z = ^ TO X.+rPY (1) 

1=1 j=i lii ^ ^ 

m 

^ ^ X 1 j=l,2j»**,n ( ^ ) 

1=1 

n 

^ X < nY 1=1,2,..., m (3) 

J=1 ^ 

X > 0 Vi and j (4) 

• X . J 

Y^ = (0, 1) V 1 (5) 

Let Kq, and Z 2 represent the sets of warehouses which 
correspond to fixed close, fixed open and free warehouse, 
then we have , 

Kq = [is = 0] 

= [1; = 1] 

K 2 = [is T 3 _ unassigned] 
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and , 

Kq UK^ UK^ = [1,2,..., n] 

The linear programming problem can now be stated as 


Minimize 

m 

z = ^ 

1=1 

n 

1 C 
0=1 

m 

X 4 . ^ p Y 

1=1 ^ ^ 

(6) 

S. t. 

in 

= 1 

t] ■ — ly^j***^ ][X 

(7) 


n 

j=i 

0 

1 c 

(8) 


n 

J=1 ^ 

n 

1 e 

(9) 


n 

r ^ 

a=i 

nY^ 

1 e 

(10) 


X > 0 
ij - 


for V' 1 and j 

(11) 




1 £ 

(12) 


Since 0, from (8) we get = 0 for all j belonging 

to the set K^. 

Furthermore, using (?)> we obtain, 
n m n 

I Z 2 = Z 1 = n 

0=1 1=1 0 = 1 = 
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Since 


n 

£ 2: 
3=1 


13 


m n 

< r r X 

~ iti 3=1 


n 

fl = n 
3=1 


the constraint (9) is redundant. 

Referring to expression (6), the first term, 


m n 

x“ r 0 X 
13 


1=1 3=1 


3 13 


^ n 

I rex 

ieX^UX^ 3"^ 


Since . 


2 =0 for 3 e E( 

-1- J 


0 


and 


m 


h = 


1=1 


le 




X 


13 


Since = 0 for 1 e Kq, the second term of expression 
(6) can be written as, 


m 


/I hh = h I 


1=1 


leX^ 


xeK^ 


How the problem given by (6) to (12) can be stated ass 



71 


Minimize Z = 




n 


S.t. 


r 0, 2 + ) p V (13 

leK^UlL^ 0=1 leK^JIv^ ^ ^ 


Z. 1 = 1 1=1,2,.,., m (14) 

ieK3_UK2 ^ 


n 


oil 


1 e K, 


(15) 


\o ^ " 


> 0 


1 £ UKg 

0 ^ 1 y Z , . . . , n 

1 £ K„ 


(16) 

(17) 


where 


A 


z = z + L F 

leK^ ^ 


The above stated problem can be further simplified by 
using the following axioms. 

Axiom 1: 


let and Y'*' represent the minimum feasible 

10 1 

solution for the problem given in (15) to (17) and if all 


> 0 then, 


n 

^ A 


0=1 


10 


= nY. 


1 £ E, 


(18) 


Yi Si 


1 £ K, 


(19) 
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Proof i 


To prove (18), notice if for some i e K, 


<1 nY , then by letting I equal to 


— r T "the value of Z can be reduced, 
n 1 j 

> 

Conseqixently for and Y^ to be optimum (18) 
must hold. To prove (19) we observe that, 


1 ^ 1 ^ 
V T' -v/- -i- r~ 


1 1 

< ~ T '£ = - r 1 = ^ 

2=1 ^ J=1 1£%UK2 ^ d=l 


Applying Axiom 1, problem (13) to 17) can be reduced to 


Minimize Z = 


-‘■J- 

3=1 leK^UK^ 


C Z. + > 

^3 ife. 


P Y 
1 1 


i£K^UZ2 


1 3 = 1, 2, . . . , n (21) 


^ ho = “ ^ ® h 

J 


( 22 ) 


Z > 0 1 e Z.UKp 

J- J 


3 — 1,2,».«, n 


Y^ > 0 1 e Z 2 


(24) 
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A further simplification of the prohlem can be 
achieved by making use of following axiom. 

Axiom 2 ; 

Define a new term g where 

^1 

^ 0 if 1 £ r 

Given any feasible solution X Y to problem represented 

J_ J X 

by (20) to (24), then 


r fY = ? r fix 

ieK2 ^ ^ 0=1 1£E;^UK2 


Proof ; 


Let X . and Y be a feasible solution to the problem 

J. J X 

given by (20) to ( ?4) , then from (22), 


I f Y = r 
1£K2 ieK2 


= L 


ieK^UK2 


y 


1 

n 

i 

n 

L 

!i 

n 

T 

n 


n 

r 


c. 

.1=1 

n 


X ) 

-IT' 


X 


Id 
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Thus the problem 
Minimise Z = 


can 

n 

'i 


3=1 


now be stated as. 


isk.ue;. 




-fL 

13 


( 25 ) 


S.t. Z 2 = 1 , j = 1,2 n (26) 

1£E^UZ2 

>0 IE JJK^ (27) 

3 — 1,2, •, n ( 

It IS interesting to note that problem given in (25) to (27) 
IS essentially an assignment problem. Thus, we have 



1 


If C 


13 



Min 

KeK^UK2 




= 0 


Jrom the definition of Eq end E-j^, and using expression 
( 22 ), we have, 

= 0 for 1 e £q 

n 

= y X /n for i e Ep 
3=1 

= 1 for 1 e E^ 



APPENDIX II 


INPOMATION POR USERS' OP TiLS PROGRAM 


This program implements the hranch and hound algorithm 
for solving the uncapacited warehouse location problem under 
probabilistic demand. The program utilizes eight heuristic 
rules for node selection, A discussion of the rules and 
their efficacy in arriving at a proper solution la given 
in the mam body of the text of this thesis. 

The data structures used are one and two dimensional 
arrays. The status of warehouses is represented by three 
integers -1, 0, 1 which stand for closed, free and open 
warehouses respectively. The program in its present form 
can be used to solve problems involving a maximum of 50 
warehouses and 100 customers. However, every care has been 
taken to make the program comprehensive. The program is 
execution time dimensioned and can be used to solve a prob- 
lem of any size depending on core memory availability simply 
by changing the dimension statements in the mam program. 



76 


Input 

Card 

Preparati 

J7o. 

on; 

Eo rmat 

Symtol 

Description 

1 


^2 

xTOPOB 

liTumber of problems 

2 


13,14 

M, N 

M ; lo. of warehouses 

N ; No, of customers 

3 


10 F. 4.2 

F( J) 

Minimum level of confidenc< 

with which j-th customer's 

actual demand must be 

satisfied 

4 


10 F 7.3 

MU( J) 

Mean of random variable 

E( 1 ) 

5 


10 F 7.3 

SI&MA( J) 

Standard deviation of the 

random variable R( J) 

6 


10 F 6.2 

F(I) 

Fired charge resulting 

from establishing a 

warehouse at site I 

7 


10 E 7.1 

0 ( 1 , J) 

Unit transportation cost 

from I-th warehouse to if 


J-th customer. Vhen 
route is prohi'bited put 
G(I,J) = l.OE 30 IS same 


format 
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Output ; 

JJ'or each problem the program outputs the following 
quantitic s s 

1. The status of warehouses (i.e. whether open, close 
or free) for each orancning rule. 

2. The set of customers supplied by each warehouses. 

3. Value of the objectiTe function for each branching 
rule. 

4. Execution time for the application of each branching 
rule (tins IS done by calling library subroutine TIME). 

5. The minimum value of objective function from the set 
of value in (3) above. 

6. The rules which give minimum objective function 
value . 

The above is a minimum output configuration. If a 
detailed knowledge of the status of the warehouses after 
each simplification is required or if the changes in the 
cost matrix are needed (This being the main working area) 
print statements can he inserted at the appropriate points. 



APPENDIX III 


PEOGEM 


LISTING 



ePANCH ANC BCUNC FRCCECUPE FCP V^ApEFCUSE LOCATICN PROBLEM 


SiejcB 
$IBFTC MAIN 
C 
C 
C 

C M = M«BEP CF PCSSIBLE kAREFOLSE SITES 

C 

C ^=TCTAL ^L^EER CF CLSTCr EPS CC EMANC FGIMTS) 

C 

C R{J)=DEMAhC CF CISTCMEP J A NORMALLY CISTPIBUTEC RANDOM VARIABLE 

C 

C AML{J)=MEAN CF TF E RANCCM VARIABLE P(J) 

C 

C SIGMA CJ5 = ST ANCAPC CEVIATICN CF THE PA|^CcM VARIABLE RIJ) 

C 

C Xtl ,J)=FPACTICN CF THE CEMANC OF CUSTCMER I WHICH Ig SATISFIED BY 

C TFE VAR CLSE LCGATEC AT SITE J, 1=1, ,M, J = l,— — ,N. 

C 

C F{I)=FIXEC charge RESLLTING FROM ESTABLISHING A WAREHOLSE AT 1 

C 

C CfI,J)=CCsT CF SUPPLING THE ENTIRE CEMANC OF CUSTOMER J FROM A 

C VAREHOUSE ICCATEC AT SITE I 

C 

C V(I)=1 IF A VAFEHfUSE IS LCCATEC AT SITE I 0 OTHERWISE 

C 

C ■ NI(n=THE SET CF WAREHCUSES THAT CAN SUPPLY CUSTOMER J 
C 

C P{J)=MINIMIM LEVEL CF CCNFICENCE WITH WHICH J TH. CUSTCMER S 

C ACTUAL DEMANC MUST BE SAISFIEC 

C 

C D(J)=DEMANC CF GLSTCMEP J MUCH ACTUALLY SATISFIED 

C 

C G!I) = F(I) WHEN NKUJIS EGUAL TO 0,0 IF KKI 1) IS EQUAL TC 1 

C 

C XFt J)=ABSCISSA ASSCCIATEC WITH THE LEFT TAiL OF THE NORMAL 

C PRgEABILTY CENSITY function FCR the CCNFICENCE level PCJ) 

C 

c 

G 

CCMMCN NCCE ,HC,CPEN 

CIMENSICN ACI50 ), AMUnCO J,CC50,100 ),CC(5C,1C0 J , C ( TOO ) tOELC 50 ) , 
2C!“EGAt5C) ,FFIf5C) ,FHJ{50),ZTEMPC8) 

LCGICAL NCCE, OPEN 
REAL LARGE 

. ncde=.trle. 

HC=I.0E3C 

GXXXXXX)<X) XXXXXXX XXXXXXXXXXXXXXXX>XXXXXXXXXXXXXXXXXXXXXXXXX5<X>XXXXXXXXX 

C ^ ' ■ , , ^ 

G I N PIT C A T A 

c ; : 

C>XXXXX>XXXXXX>XX>XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 



BRANCH ANC BCUNC FRCCECURE FCR WAREHOUSE ICCATICN PROELEN 


REAC(5,1CC3NFPCE 
ICO FCRNfiT (12) 

DC 11 NPPCE=1,NFPCE 
REAi:t5,lCl)N,N 
ICl FCPNAT(I3,IA) 

READ (5 ,1C2) (PC J), J=1 ,N} 

1C2 FCRMT{1CFA.2) 

REACC5,1C3) 

1C3 FCRNATdCFl .3) 

READ f5,1C4) {SIGNftI 
1C4 FCRKATdCFT.B) 

REAClSaCS) (FU )»I = ltF) 

1C5 FCRNAT(1CF^.2) 

REACt5,1C6HfC{I,J) 

1C6 FCRNATdCET.l) 

C >X)<)(X)(>)0)<X>f X>>X>>)(X>)<X>>t)(XXXX><X>XXXXXXXXXXXXXXXXXXXXXXX)()iXX)<XX)()<X>X>X 
CALL CCS! {F,C,C,P,APU,SIC-f'A»XP,CC,P,N,KK ) 

CALL TINE d PL) 

NRLLE=1 

IF {NPLLE.EC .1) GC TC 5CC1 
FOCI feRnE{6,5GC9) 

5CC9 FCRNAT {/ ,ACX,=fPPCBLEN IS BEING SOLVED EY LARGEST DELTA RLLE=*) 
),RITE 16,5011) 

5017 FCRNAT (35X ,56aF-) ) 

22 CALL CTSIP (KK,F,CNEGA,AC,DEL,C,M,N1 ) 

IF(NRLLE.LE.4) CF£N=.TRUE. 

IF {NRULE.GT.4) CFEN=.FALSE- 
DC 23 1 = 1, N 

IFdKtl ) .EC.O) GC TC 24 

23 CCNTINUE 
GC TC 70 

24 II=C 

CALL NRULES I CN EG A , CEL , C , Y , X t G ,CC ,C , F, KK, PHI , NRULE ,H,N, PH J ) 

DC 3C 1 = 1, N 

IFIXKtn.NE.O) C-C TC 30 
11=1 ■ 

3C CCNTINLE 

IFfll.EC.O) GC TC 70 
NCCE=. FALSE. 

GC TC 22 

70 CALL PRLIN (6 , X ,Y ,cC ,C,F tHH , h, Z ) 

VvPIT£t6 ,3CC1>Z 

3Gei FCRFATf / ,15X ,=4*VALUE CF CEJEC. IVE FUNCTICN = *,E15 .8 ) 

ZTENF { NPLLE ) = Z 
DC £ 1 = ] ,N 

4CC1 FCRf“ATt/,9X,I2,5X,8F12.1) 

8 CCNTINLE 

CALI TIRE INPL) 

NPLJLE=NRLLE4l 

IFCNRULe.EG.2) )^R ITE 1 6 ,5010 ) 

IFCNRULE.EC.3) R R ITE (6 ,5C1 1 ) 



o o 


BPahCH iSKC OCUNC PRCCECLPE FCP WAREHOUSE LOC/iTICN PROetEf' 


IFfhRULE.EC .4) ^PITEf6,5012> 

IF CNRLae.EC .5) i^RITE 16,5013) 

IF.^^RULE.E£ .6) W P ITE ( 6 , 501 4 ) 

IF{^RULE.EC.7) WR ITE 1 6 , 50 1 5 ) 

I FfhPULE.EC .8) )sRITE{6,5Ql6) 

ECIC FCRN4T {/ ,4C>,«FFC8LEM IS EEING SOLVED EV L/SRGEST Of'EcA RULE*) 
soil FCRRiST(/,4CX,*FPCBLER IS EEING SOLVED BY LARGEST DEMAND RULE*) 
5C12 FCRNiiT(/,4CX,*PPCBLER IS BEING SOLVED BY LARGEST Y RULE*) 

5013 FCRNAT (/ ,4C>,*PPC8LEf' IS BEING SOLVED BY SMALLEST DELTA RULE*) 

5014 FCRNAT(/,4GX,*PFCeLEN IS BEING SOLVED BY SMALLEST CMEGA RULE*) 

5015 FOR^ATf / ,4CX,*PFCBLEM IS EEING SOLVED BY SMALLEST DEMAND RULE*) 

5016 FCRNAT(/,4C>,*PPCEIEM IS BEING SOLVED BY SMALLEST Y RULE*) 

kRnE{6,5cn) 

NCCE=.TRIE. 

DO 15 I=1,N 
KK(I)=C 
CO 15 J=1,N 
7^5 ca ,J)=CC(I ,J) 

' IF(NPULe.GT.8) go to 80 
GO 10 22 

60 CALI AMINCE tZTENF) 

CALL TINE (NTPC) 

11 CCMINIE 
STOP 
END 

fiBFTC COST 

C ******** ********’!***************************** *****#3>3{t******' *#*!*<****: **5 

SLBPCUTINE COST FOR CALCULATING Cl I ,J ) MATRIX 

SLBRCUTINE COST 1 F , C , C , F , AMU, S IGM A, XP, CO, M , N, KK ) 

DIMENSION C (M,N),C(N) ,F IN ) , AMU ( N ) ,S IGM A IN ) ,XP IN ) , F ( M ) ,CG(M,N ) , 
IKHIM) 

DOUBLE PFECiSiCN Z2,CISTT 
CCDE=1.0 
DC 1C v]=l,N 
p I S 1 T = p I 4 ) 

CALL NORMAL IZZfCISTT, CODE) 

XFt J)=ZZ 

DTJ)=SIGMAlu)*XF{J)+AMUIJ) 

1C CCMINLE 
DC 6 1 = 1 ,M 
6 CCNTINUE 
DC 11J = 1,N 
DC II 1 = 1 ,M 
C?I ,J) =C {I , J)*C IJ) 

TV CCNTINLE 
DC V 1=1 ,M 
5 KI<n)=C 

KPITEI6,2C13) 


BRAKCH AKC 8CUNC FRCCECliRE FCP WAFEHCUSE LOCATICN PrOBLER 


2C13 r-CRNaT{/,56>,>«'STATUS CFWAPEPCUSES*/ t56X ,1 S C 1H-) ) 
kPITE{6 t2C15) (KHd ) = 

2015 fCPPAT (/ ,49X,5I€) 

DC 6CQl=liP 
DC 6CC J=1,A 
6CO CCCI ,J5 = Cd ,J) 

RETIPK 

EPO 

$IBF7C ^CR^'AL 

C ELERClTI^E NCRPAL FCP CAlCLLATING THE VALUES OF ABSCISSA OF NORMAL 

C CISTRIBLTICN FCR GIVEN VALUES OF CONFIDENCE LEVEL 

4 : ^ # $ $ 44 # 4 # 44 44 

SLBRCUTINE NCRM AL C Z Z t C ISTT » CC CE ) 

CCLELE PPECISICN Z Z » C ISTT ,Z ,C 1ST, AL FHA 

T=C.C1 

A=C,C 

Z=C.CDGO 

CIST=C1.5CCC 

C44444:cCcE=1.C NEANS that DISTRIBUTION IS GIVEN AND z TO BE CCmPUTEO 
IFICCDE.EG. 1.0) GO TC 7 

C44444paRT-I GIVEN Z.CISTRIEUTICN TO BE FOUND 
IF{ 2Z.EG.C.C) GC TC 100 

10 IFt (ZZ-Z).GE.O.CI) GC TC 1 
T=ZZ-Z 

1 ALPHA={T/12.0)’dCF(Z)+4.0*CFCZ+0.25=»T) + 2.C*DFfZ+C.5*T) + A.0*DF(Z40. 
175=5'T) + CF (Z+T) ) 

Z=Z + T 

OIST=DIST+AlPHA 
IFCZZ.EC.Z) GC TC ICO 
GC TC 10 
ICG DISTT=C1ST 
GC TC 555 

C44444PART-II CISTRIEITICN GIVEN , Z TC BE FCLND 
7 IFICISTT.LT.O.S) GC TC 77 
IFCCISTT.EC .0.5) GC TC 1000 

11 ALPHA={T/12 .0)=>{CF( Z)+4.0«CF{Z+0.25=fT) + 2.C4CF(Z+0.5*T)+4.Q*DFCZ+0. 
175*T) + DF 1Z + T) ) 

■ 1 = 1*1 

CIST=DIST+A LFHA 
IF(CISTT.LT.DIST) GC TC ICl 
IFUISTT.EC.DIST) GC TC 102 
GC TC 11 

C444444C0NPARI NG T WITH ACCURACY CESIREC 
ICl IF(T.EC.C.CCOl) GC TC 102 
Z = Z-T 

CIST=DIST-AIPHA 
C44444SFECIFY ACCLRACY 
T=G.C001 
GC TC 11 

17 DISTT=CISTT+0.5C00 



BRAKCH BCl^C FRCCEGUPE FCR kAFEhCUSE L0C4TICN PROELE?^ 


A = 1.0 
GC 7C 11 
1C2 Z 2= Z 

IF{ «.EC.1.C ) ZZ=-ZZ 
IFCA.EC.l.C ) CISTT = CISTT-Q.5COO 
GC 7C 555 
ICCC ZZ=C.O 
555 Ret CRN 
END 

^IBFTC CF 

FINCTICN CF IZ) 

CE^SITY FUNCTICN 
CCLBLEPREC I SIcN Z 

DF=CEXF {- { Z =!*2/2.0CC0) )/CSGFTUA.ODCO/7.CCCG) 

RETIPK 

END 

$I8FTC CTSIP 
C 

C SteFCUTINE CTISN FCR c'^CLe ThPCGF SIMPLIFICATION CRITERIA AT A NODE 

C 

C ^ * :{! * ^ 4 4 * 3{ :{ ^ 5{ 3f :! 5} )|. 3}: 4 4 sj! < :5 ^ ^ !« # * ♦ * * 4! * * 5> 3}, ^ * $ « * « # # :{e « 

SlBPCyTINE CTSIN t KK » F,Cp EG A, AC, C EL ,C, N,N ) 

CCNNCN NCCE ,HC,CFEN 
LCG ICAL NCCE ,CFEN 

CINENSICN KH(P) ,F(|«) tC.PEGA(N),AC{M)tCELIM)»C(M,N) 

2013 FCRRAT {/,56X,*STATLS CFW AREFCUSES*/ ,56X, 1SC1H-) ) 

2015 FCRNAK/ ,A5X,5I6) 

IF (NCCE ) GC TC 1 
IF tCPEN ) GC TC 2 
Gc TC 3 

: 

C CRITERION AFFLIEC FCR FIX OPENING OF A NAREFOCSE 

C : — — — 

1 DC Id 1 = 1 ,N 

IFCNKdI.NE.OI GC TC 10 

SLN=C.C 

CC 20 J=lfN 

DC 30 K=1 ,N 

I FIN. EC. I) GC TC 22 

IFCHKtl) .EC.t-D) GC TC 22 

IF(C IK,J) .GE.1,GE25) GCTC22 

ACIK)=C(K,J)-C(I,J) 

IF(AC(K) .LT.0.0) GC TC 41 
GC TC 30 
22 AC(i<) = l.CE2 5 
20 CCNTINLE 

i'N=N-i : 

AMN=AC( 1 ) 

DC 35 K=1,FN ’ 

IFIAC(K+1).LT.AMN) qC TC 28 



o o o 


EPA^CH AKC ecu^c FRCCECLRE FCR WAPEI-OUSE LGCA7ICN PROBLEf* 


GC TC 35 
2£ AMF = ACf!< + l) 

3 5 CCMIME 
GC 7C 23 
41 AMh=0.0 
23 Slf' = StK+AM> 

20 CCMIME 

CEL {I) = StN-F(n 
kRnFt6,2C2C)I ,CEL (I J 
2C2C FCRFAT(/,5G>,I3,F2C.3> 

1C CCMIME 

DC AC 1=1 ,f 

IF(!<Kn ).N£.0) GC TC 40 
IFCCEMI I.Gl.C.C) GCTC 21 
AO CCMIME 
RETIPN 

21 CC 54 1 = 1 

IFtf<Kn ) .^E .0 ) C-c TC 54 
IF( CELd I.IE.O.C) GC TC 54 
KHCI)=1 
5 4 CCMIME 

VFI7E(6 ,2C13) 

^PITEI6»2C15) IKKtl) ,I=1,F) 

2 DC 55 1=1, F 

IFf FKl I) .EC .0) GC TC 60 
55 CCMIMe 
RETIPN 


CPIIERICN AfFLIEC FCR PECUCIFG TFE NO. GF NI 


60 DC 70 1=1 

IF! HKd) .FE.l) GC TC 70 
DC 80 J = 1,N 
eMF=Cd,J) 

DC GC K = 1 ,F 

IFCC(K,J).IT.ENIN) GC TC 80 
GC CCMINIE 

CC ICC i<=l,F 
IFf KKCKI .EC.I ) GC TC 100 
C1K,J)=HC 
ICO CCMINLE 
80 CCMIME 
70 CCMIME 

DC 110 1 = 1 ,F 
CC 120 J = 1 ,N 

IFfC (1, J>.n.l.0E25) GC TO 110 
120 CCNTINLE 
K}<tl) = -1 
lie CCNTIME 

DC 130 1 = 1, N 

IFlKKd ) .EC.O) GC TC 140 



BRANCH J!NC 8 CLNC FRCCECCPE FCR NAFEEOUSE LOCATICN PROELEN 


13Q CCMIME 
RETIRN 

G 

C CRlTERICN AFFLIEC FCR FIX CLCSING CF A WAREhOUSE 

C 

lAO CC 15C 1=1, F 

IF(KK{ I ) .NE .O ) GC TC 150 

SLN=C.C 

DC 1€C J = 1 ,N 

CC 17C i< = l,N 

IF(l<KtK).NE.1) GC TC 152 

IFIC {f<, J).GE.l.GE25] GC TC 152 

AC{K)=CtK,J)-C{I,J) 

IF(ACfHl.n.O.O) GC TC 162 
GC TC nc 
152 ACCH ) = 1.CE25 
17C CCMIKCE 
PP=N-1 
CMN=ACn) 

CC ISO K=1,?N 

IFIACIK + D.LT.CMM GC TC 158 
GC TC 18C 
158 CMA=AC IN + I ) 

180 CCMIME 
GC 7C 183 
lt2 CMN=C.C 
183 Slf = SUf'+CMN 
160 CCMIME 

CFEGA (I1=SIN-Ff I) 

2023 FCRNAT(/,5CX,I3,F20.3) 

150 CCMIME 

nr 138 1 = 1, IK 
IFfxKni .NE .0) GC TC 138 
IFICPEGAtD.LT.C.O) GC TC 1A3 
138 CCMIME 
RETIPN 

1A3 CC 1A4 1=1, P 

IFIKKCn.NE .0) GC TC 144 
IFtCNEGAIIJ.GE.G.O) GC TC 144 
Kl<tn=-1 
CC 167 J=1,K 

CtI,l)=HC 
167 CCMINUE 
144 CCMIME 

3 CC 2C2 1=1 rP 

iFti<K{n.M.(-rn gc tc 202 

DC 204 J = 1 ,H 
2C4 C !I ,J )=HC 
2C2 CCMIME 

VvRITECe ,20131 

NpnEt6,2C15)IKH111,I = l,P) 



BRAhCH /She BCt^C FRCCECLRE FCP KAFEFOUSE tCC/^TICN PROBLEf' 


CC 1 1=1, f 
7 CCNTIME 
CC 5C3 1 = 1 

IFfKKf n.EC.O) GC TC 201 
2C3 CCMTME 
RETIRN 

C — 

C FLR7HER S I f F I L I FI C ^ T I CK FCR REDUCING TFE NO. OF NCI) 

C — — — , 

2C1 CC 270 1 = 1, N 

IFIFKCI) .NE.1 ) C-C TC 270 
CC 28Q J=1,N 
CMN=C{I,J) 

CC 200 K=l,> 

IFCC CK ,J) .LT.CNIM GC TC 280 
200 CCMIME , 

CC ^81 K=1,F 

I FC 1<KC K) .EG .1 ) GC TC A81 
C CK ,J) =HC 
4E1 CCNTIME 
2EC CCNTIME 
210 CCNTiNUE 
GC IC 1 
END 

fciBFTC FRIIN 

^ iff ^****t** ******** ** 

C 

C SLBPCUTINE FRLIN.FCR SCIUTICN CF LINEAR PROGRAMNINC- 

C 

if i)i)i!**iji**!ii*** ****** ii*****:^** *** ************************************ ** 

SIBFCUTINE FRLIN{G,X,'r,CCtC,F,KK,f',N,Z) 

CCNNCN NCCE ,HC,CPEN 
LCGICAl NCCE ,CFEN 

DIf'ENSICN G CP) ,X(N,N),YCP),CCCM,N ) ,F C N ) , KK CM ) , C( H,N ) 

CC 15 J=1,N 
CC 16 1 = 1, N 

IFCG Cl , J ).IT.1.CE25) GC TO 15 
16 CCMIME 
ivPITECE ,7) 

7 FCRNATC2CX,=«SCLITIcN IS INFESIeLE*) 

Z=HC 
RETLRN 
15 CCNTINUE 
DC 1C 1=1, N 

TFCKKCI).EC.(-l)) GC TC 10 
IFCNK(I).EG.l) GCn=0. 

IFCNNCD.EG.O) GCn = FCn 
1C CCNTIME 

CC 20 J=1,N 

ENTN=HC 

CC 30 i< = l ,N 



eRA^CH ANC BCUKC FPCCECURE FCR WAREHOUSE LOCATION PROELEH 


IFfl-KCKI.EG.T-l)) GO TC 30 
CALL NCUST IM,K,C,H,N) 

IFC iCC(!<,J)4G{H)/FLCATIMn.G£.Ef'IN) GO TO 30 
KVAPE=K 

EMK=CC ?K ,C) + G (K5/FLCATIM) 

30 CCMINLE 

>;IKMRE,J5 = 1.0 
DC AC 1 = 1, F 

IF( I.EC.HVAPE J GC TC 40 
XM , J)=O.C 
AC CCMInLE 
2C CCMINLf 

z=c.c 

DC EC 1=1 

ifikkt n.Ec.f-in v{i)=g.o 
IFIKKI D.EC.l) T(n=1.0 
IFCHKII I.EC.OJ GC TC 55 
GC TC 65 
55 ASUT'=C.C 

DC 60 L=1 ,N 

60 ASU^=ASL^ + >< (I ,L) 

CALL NCUST 1KJ,I,C,H,M 
Yd )=ASLNyFLCAT (KJ) 

65 CC TC J=1 

1C Z = Z4CCfI,J)=*XtI,J) 

2=Z+FfI)«TlI) 

50 CCNTIME 
RETURN 
END 

$ieFTC NCIST 


C SLBPCUTINE NCLST FCp CALCUL AT ING TH E NC.OF NICU 


SLBPCUTINE NCLST UZ,HL,C,H,N3 


CCPNCN NCCE,hC,CFEN 
LCGICAL NCCE,CPEN 
CIFEfsSICNC t NtM 
■ NZ=C 

DC 10 J=1,N 

IFIC {KL,J3 .IT.HC) NZ=NZ+1 
10 CCNTIME 
RETURN 
END 

$IBFTC NRLLES 

c 

C SLBFCUTINE NRLLES FCP AFPLICATICN OF BpANCHING RULES 

c : 

4 :{1 ^ >J :{ 3{ 4 :} 4 :5 ^ 3} :f 4 jJ li :} 3} 3534 ^ 4 14:4 34^ 34:4:4:34 4 # 34 34 34 344 444 4 4444 44:34 4 434 44 44 

S LB F CUT I N E N RULES ( GF E GA , C El » C , Y , X, G, CC ,C , F ,KK, PHI , NRUL E,N , N , PH J 3 
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BRANCH 6NC BCUNC FRCCECLFE fCR ViAREHOUSE LGCATICN PROELEH 


CI^ENSICN FMIR),CNEGAtR} ,CEI(F),C|N),Y{h }»XCM,N),GIM),cC*^»N 5* 
1C (F ,N) ,F(F) jHHfM . 

INTEGER RILE 
CC 1CRLLE=NFULE ,8 

IF(FULE.EC..1.CR.RULE.EC.5) C ALL CY Cl CEL, PH I , M, N, KK, RUL E ) 
IFIRLLE.EG.Z.CR.RLLE.EC.S) CALL CYC f CREC-A , PHI ,P,N,KK, RULE J 
iFt RULE. EC. 3. CR. RULE. EG. 71 CALL CY C( C, FF I,R,N,KK,RLLE) 
IF(RLLE.EC.A.CR.RLLE.EC.81 GC TC 2C 
15 IF(PLLE.LE.A) call AFAX fPH,KK,F,N,RLLE,C,PHJ ) 

IFfFLLE.GT.A) QfiLl AFiN f PH ,KK , F, N, RUl E , C , PHJ 1 
RETURN 

20 CALI FRLIN (G , X , Y , CC ,C , F , XK ,F , N , Z 1 
CALL CYCn,FHI,N,N, XK, PULE) 

GC TC 15 
1C CCMIME 
ENC 

$ieFTG CVC 

SLB FCUTINE CYC ( FS I , FF I , F ,N , X K , RULE 1 
CIFENSICN FSMF),PFI(F),KK(F) 

INTEGER RLLE 

IF (RULE. EG. 3. CP. PULE. EC. 7) GC TO 15 
DC IC 1=1 ,N 

IF (XKd l.NE.Q) GC TCIO 
PFI (I)=PSl (j) 

10 CCNTINLE 
GC TC 20 
15 CC 18 1 = 1, F 
18 PFI {n=PSI (I) 

2C RETLPN 
END 

$I8FTC AFAX 


BRANCHING RULES FCR FIXED CFEMNG OF A kAREFOUSE 


SIBPCLTINE AFAX{FHI,XX,F,N,FULE,C,PFJ ) 

CIFENSICN FFI {F),XK(F),C{F,N),PFJIM) 

INTEGER. RULE 

IF(RULE.NE.2) GC TC 7 

CC 6 I=1,F 

PFJ (I)=O.C 

IFtXXfll.NE.C) GC TC 6 
DC 5 J=1,N 

IFCC(I,J1.GT.1.0E25) GC TC 5 
PHJ(I1=PFJ(1) + PFHJ) 

5 CCNTINLE 

6 CCNTINLE 

CC 12 1=1, N 
12 PFI 111= Pf 4(1) 

7 RFA>=-1.CE26 
II = C 

DC 1C 1=1, N 



o r> r> 


BRAhCH AKC BCUNC FRCCECIRE FCR V^APEHOUSE LOCATICN PROBLEM 


IF tKKlI 3 .hE.O) GC TC 10 
IFCFFI (D.GE.PMAX) GC TC 15 
GC TC 1C 
15 RRA>=PHI ( I) 

11 = I 

IG CCKTIKLE 

IFdI.EG-C) ^PITEC6,2C) 

HHt 113=1 
RETIRN 

2C FCRPAT (IFC/Z’^AL L VALUES LESS THAN -1.0e36*1 
END 

$IEFTC AM IK. 

C— — 

C BRAKCHIKG PILES FCP FiXEC CLCSIKG OF A WAREHOUSE 

— 

SLBRCLTIKE AMIK C FH I , K K , M ,K ,PUl E, C» PH J 3 

CIMEKSICK FHI(M5,KK(M3 ,C IN ,K )» FHJt N3 

IKTEGER RULE 

IF{PLLE.KE.13 GC TC 7 

CC i I=l,t 

FHJII3=n.C 

rFIKKCn.KE.OI GC TC 6 
• CCS J=1,K 

IFICdf J3.GT.1.GE25) GC TC 5 
PH J (II=FH Jt I3 + FHH JJ 

5 CCKTINUE 

6 CCKTIKlE 

CC 12 1 = 1, K 
12 PHI CI3=PHv 1 13 

7 RMK = 1.0E3€ 

HC=1.0E3C 
II = C 

DC 1C 1 = 1, K 

IFIFKin .KE.Q 3 GC TC 10 
IFIFHI {I5.1E.RMIK3 GC TC 15 
GC TC 10 
15 RMIK=PHI{I3 
11= I 

10 CCKTIKUE 

IF{1I.EQ.C3 PRINT 20 
TF{ II.EG.C3 WRITEt6,2C) 

KK{II3=-1 
DC 25 J=1,K 
25 CnT,J3=HC 

2C FCRMAT(1HC//*ALL VALUES GT 1 . 0 E 26 «) 

RETURN 

END 

le FTC AM IKCB 


SLBPCUTIKE AMINCE FCR FIKCIKG MINIMUM CEJECTIVE FUNCTION VALUES 



{/I iti ru 
LTt o Kn 


eRA^CH A^C BCUNC FRCCECLRE FCP WAFEhCUSE LOCATICN PROBLEM 


SlBFCUTI^E iNIKCR (ZTERF) 

CI^'E^SIC^ ZTEKF (8) ,ITEPPC8) 

RMN = 1.CE36 
DC 1C L=l,8 

IFf ZTERP tL) .GT.FRU} GC TO 9 
RM h = ZTERF { L) 
ne^F^L) = L 
GC TC 10 
9 nE^F(t)=C 
1C CCMIME 

VvPnE ( 6,20 5 PRIK 

2C FCRR^Tdh// ,20X VALLE OF OBJECTIVE FUNCTION IS ==«,2X,ei5. 

13) 

ViRnE{6,55 ) 

CC5C 1=1,8 

IFCITERPd) .EG.C) GC TC 50 
feRITH ( 6 ,25 5 ITERF(L ) 

FCBRATdF// ,20X,I2) 

CCMINLE 

FCRRSTt/ ,2CX,=»TFE FCLLCV^ING RULES GIVES TFE MlNiNUR OBJECTIVE FUNCTI 
ICR*) 

RETLRN 

ERD' 

$ENTRY 



